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The classical spin-flop is the field-driven first-order reorientation transition in easy-axis anti- 
ferromagnets. A comprehensive phenomenological theory of easy-axis antiferromagnets displaying 
spin-flops is developed. It is shown how the hierarchy of magnetic coupling strengths in these anti- 
ferromagnets causes a strongly pronounced two-scale character in their magnetic phase structure. In 
contrast to the major part of the magnetic phase diagram, these antiferromagnets near the spin-flop 
region are described by an effective model akin to uniaxial ferromagnets. For a consistent theoreti- 
cal description both higher-order anisotropy contributions and dipolar stray-fields have to be taken 
into account near the spin-flop. In particular, thermodynamically stable multidomain states exist 
in the spin-flop region, owing to the phase coexistence at this first-order transition. For this region, 
equilibrium spin-configurations and parameters of the multidomain states are derived as functions of 
the external magnetic field. The components of the magnetic susceptibility tensor are calculated for 
homogeneous and multidomain states in the vicinity of the spin-flop. The remarkable anomalies in 
these measurable quantities provide an efficient method to investigate magnetic states and to deter- 
mine materials parameters in bulk and confined antiferromagnets, as well as in nanoscale synthetic 
antiferromagnets. The method is demonstrated for experimental data on the magnetic properties 
near the spin-flop region in the orthorhombic layered antiferromagnet (C2H5NH3)2CuCl4. 
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I. INTRODUCTION 

In antiferromagnetic crystals with a preferable direc- 
tion of the magnetization, a sufficiently strong mag- 
netic field applied along this easy axis direction "over- 
turns" the sublattice magnetization vectors Mi and M2 
(Fig. 2J. Neel demonstrated this threshold-field ef- 
fect theoretically for classical two-sublattice antiferro- 
magnets with sufficiently weak anisotropy in 1936i Fif- 
teen years later, this prediction of a jump-like reorien- 
tation transition driven by an external magnetic field 
was confirmed by experiments on the antiferromagnet 
CuCl2-2H20. 2 Since that time the spin-flop transition 3 
has been observed and was investigated in great de- 
tail for a large group of antiferromagnets (see, e.g., 
Refs. l ^'lGllTip a nd further examples and bibliography 
in Refs. |9ll0lll| ). Originally, the name spin-flop (SF) 
transition was restricted to the field-driven reorienta- 
tion transition in two-sublattice collinear easy-axis an- 
tiferromagnets following Neel's prediction. However, in 
many other classes of antiferromagnets similar mecha- 
nisms cause different types of field-driven reorientation 
effects so that one can speak about a class of spin-flop 
phenomena in a wider sense. These include field-driven 
transitions in antiferromagnets with a Dzyaloshinskii- 
Moriya interactionj&i. 2 in multisublatticeji 3 ^ and in 
noncentrosymmetric antiferromagnetsii 

The spin-flop transition comprises in a simple form 
main features of magnetic reorientation and phase tran- 
sitions and it gives rise to various remarkable physi- 
cal anomalies. It was found that the spin-flop transi- 
tion is very sensitive to direction and values of the ap- 



plied field 15-16 and specific transitional domain struc- 
tures are formed in the spin-flop regioniiSiii These ef- 
fects have been observed in MnF^A^ and some other 
antiferromagnetsiiSiS . Theory also predicts a specific re- 
construction of the domain walls in the vicinity of the 
spin-flopi2i2i and extraordinary surface effects, which 
have been called "surface spin-flop" M> The investigations 
on spin-flops made important contributions not only to 
different research topics in magnetism, but also to such 
general fields of physics as thermodynamicsjiS*^ 3 . nonlin- 
ear physics ^^iS 5 . and the theory of phase transitions and 
critical phenomenal 
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FIG. 1: Basic spin-configurations in collinear uniaxial an- 
tiferromagnets. At zero field the magnetization vectors mi 
are antiparallel (a) . The magnetic fields along the easy direc- 
tion H > Ho stabilizes the spin-flop phase (b). In magnetic 
fields deviating from the easy direction canted states are re- 
alized (c). In a sufficiently strong magnetic fields, H > H c , 
the spin-flop and canted states transform into the saturated 
flip phase (d). These solutions are degenerate with respect 
to the sign of the staggered vector 1, which is represented by 
two-headed arrows in the bottom panel. 

During the last decade spin-flops have been iden- 
tified and studied in many novel bulk antiferro- 
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magnets including such classes as antiferromagnetic 
semiconductors j2i organic magnets^ in cuprates, such 
as La2CuC>4 and Nd2CuC>4, being base materials for 
high-temperature superconductors^ 3 -, or in noncen- 
trosymmetric antiferromagnets^iSS However, the focus 
of current interest in spin-flops has now shifted to- 
wards nanomagnetic systems. Many of recently syn- 
thesized nanostructured materials have magnetic con- 
stituents with antiferromagnetic coupling. This vast 
and intensively investigated class of antiferromagnetic 
nanostructures includes different types of ferromag- 
netic/antiferromagnetic bilayers^ Synthetic or artificial 
antifcrromagnets are recently designed for high-density 
storage technologies, spin-valves and magnetic random 
access memory (MRAM) devices i^i* 3 ^ In addition to lay- 
ered systems nanoparticles of antiferromagnetic materials 
are currently investigated, 3 ' A variety of size- dependent 
electronic and magnetic effects are intensively investi- 
gated in many of these nanostructures ,24^ In particu- 
lar, experiments indicate that in ferro/antiferromagnetic 
bilayers the reorientation within the antiferromagnetic 
structures strongly influences the interface interactions 
and thus magnetic properties of the ferromagnetic 
subsystem. ] 36 ! 37 A rich variety of specific field-induced re- 
orientational effects have been found in synthetic antifcr- 
romagnets such as toggle mode MRAM devices^S* 3 ^ and 
finite antiferromagnetic superlatticesi 3 ^*^ 

In a broader context, new methods for the investi- 
gations of antifcrromagnets in the bulk or close to an- 
tiferromagnetic surfaces^idii^ and of nanostructured 
antiferromagnets 3 ^ 3 ^ 6 ^^ provide a wealth of data on 
spin-reorientation effects. These investigations require a 
refined theoretical description of these phenomena. The 
theoretical understanding of the the spin-flop was de- 
veloped during many years of investigations with con- 
tributions from many researchers and great progress to- 
wards an understanding of the magnetic states near the 
spin-flop transitions has been reached. Still, the present 
stage of the theory is mostly dedicated to separate phe- 
nomena and a coherent description of the magnetization 
processes near spin-flops was not presented yet. This 
problem is evidently one that should be addressed with 
phenomenological approaches of magnetism, i.e., a "mi- 
cromagnetism" for antiferromagnetics. 43 First solutions 
were obtained for fields applied along the easy axis^ 4 ^. 
The stability limits forming "critical astroids" in oblique 
magnetic fields, i.e. in fields deviating from the easy-axis, 
and the "critical angle of the spin-flop" were calculated 
in Refs. 0] for a model with second-order anistropy. 
The influence of the fourth-order anisotropy on the 
metastable states has been investigated in Ref. j^. In 
Ref. 0| the problem has been reduced to an effective 
model akin to the usual reorientation in a uniaxial ferro- 
magnet, and a complete set of solutions was derived. 

In this paper we give a comprehensive phenomenolog- 
ical theory of the magnetic states and their evolution in 
applied fields for a two-sublattice collinear antiferromag- 
net. The model is explained in section [HI We calculate 



all possible magnetic configurations in the system (sec- 
tion IIII(I and give a physically clear description of the 
main features of the magnetization processes. In partic- 
ular we calculate and analyse the tensor of the static mag- 
netic susceptibility (sections IIII Bl and IIV|I . This enable 
us to generalize and systematize results on bulk spin-flop 
and formulate directions for the investigation of this tran- 
sition in bulk antifcrromagnets and in antiferromagnetic 
nanostructures. In section IIVI the occurrence of mul- 
tidomain states by demagnetization effects is analysed. 
There are characteristic peculiarities of the field and an- 
gular dependencies of the magnetic susceptibility, which 
can be employed in experiments on the spin-flop transi- 
tion. In sectionEJwe demonstrate this approach based on 
experimental data for orthorhombic (^HsNHa^CuCLi, 
which is a layered model antiferromagnet^ 7 - 

II. MODEL AND EQUATIONS 

Within the phenomenological theory of magnetism 
the magnetic (free) energy for a bulk collinear two- 
sublattices antiferromagnet can be written in the follow- 
ing formifl*^ 

W = J to(mi,m 2 )dV (1) 

= J | Jmi • m 2 + e (mi,m 2 ) 
-H(°) • (mj + m 2 )M- 1 
-iH< m ' • ( mi + m 2 )M^M^dV. 

We assume here that the vectors of the sublattice magne- 
tizations Mj do not change their modulus and their ori- 
entations are described by unity vectors uij = Mj/Mo, 
Mo = |Mj|. Hence, we develop our theory for constant 
low temperatures in the antiferromagnetically ordered 
state. The energy consists of the exchange inter- 
action with exchange constant J, the magnetocrystalline 
anisotropy energy e a , and Zeeman energy contributions 
due to the external magnetic field H^ c ) and the demag- 
netization field H", the latter giving the dipolar stray 
field energy. In this paper we are interested in anti- 
fcrromagnets with a preferable direction of a magnetic 
ordering, i.e. easy-axis systems. In uniaxial antifcrro- 
magnets the easy-axis coincides with the principal axis 
of symmetry, and in orthorhombic crystals with one of 
the orthorhombic axis. In this paper we chose coordi- 
nates such that the easy axis is along the z-axis. The 
energy density e a — ea + ei includes uniaxial e« and 
in-plane e£ parts. The uniaxial anisotropy e£ consists 
of invariants of type Tn±™ where n is an integer, 

and k = l,2..nM These energy contributions strongly 
decrease with increasing n. In the vicinity of the spin- 
flop field the second-order uniaxial anisotropy is "can- 
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celled" by the applied magnetic field. Hence, higher- 
order terms must be be taken into account 1^2^ There- 
fore, both second-order and fourth-order terms must be 
included in the theory for spin- flops. One can write the 
uniaxial anisotropy asifi 



e { a ) {m lzi m 2z ) = 



K 

-^( m lz + m iz) ~ K'm lz m 2z 



K20 
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K21 



(mfa + m 2a) 

m la m 2z (m 2 la + m 2 a ) 



K22 2 2 

— m la m 2a- 



(2) 



Usually, the second-order terms with constants K, K 1 
play the dominant role for the orientation of the magnetic 
vectors. But, the fourth-order terms with n = 2 and 
constants K 2 k become vital near the spin-flop, as will be 
shown later. 

The anisotropy energy in the basal plane includes 
invariants composed of in-plane components of the vec- 
tors m.i and, depending on the crystal symmetry, includes 
quadratic terms for orthorhombic crystals or terms of 
higher degree for uniaxial crystals. These energy con- 
tributions lead to additional orientational effects. They 
determine one preferable direction for orthorhombic sym- 
metry or more directions for uniaxial antiferromagnets 
for the spin orientation in the basal plane perpendicular 
to the magnetic easy axis. 

The model Eq. Q with anisotropy @ describes a 
vast group of antiferromagnetic crystals with collincar 
order, including such well-studied compounds as MnF2, 
FeF 2 , Cr 2 3 , GdA10 3 and others^ In this class of an- 
tiferromagnets effects of magnetic couplings are absent 
that violate a collinear and spatially homogeneous or- 
der in the ground state, such as competing exchange or 
Dzyaloshinkii-Moriya interactionsiii 

The equations minimizing the energy functional 
together with the Maxwell equations for the magneto- 
static problem determine the distributions of the mag- 
netization fields rxii(r) and the stray field H^ m '(r) in 
the sampled These integro-differential equations are too 
complex and impractical, even for brute-force numerical 
calculations. However, the problem can be reduced to 
a number of simplified auxiliary problems Follow- 
ing these standard methods we start from the analysis of 
spatially homogeneous states in fixed internal magnetic 
fields H = H( c ' +H( m ). 

In the antiferromagnetic crystals displaying a spin- 
flop, the exchange coupling is much stronger than the 
anisotropy energy, J 3> w a . In this case it is convenient 
to use the net magnetization vector, m = (mi + ni2)/2, 
and the staggered magnetization vector (or vector of an- 
tiferromagnetic order), 1 = (mi — rri2)/2, as internal vari- 
ables of the system. Because |nii| = 1 these vectors sat- 
isfy the constraints m • 1 = 0, m 2 + l 2 = l4i Indepen- 
dent minimization of Q with respect to m yields (see 



Refs. ITlllT^ for details) 

m = [H — n(H n)]/JT , H c = (2JM ) , (3) 

where n = 1/|1| is the unity vector along the staggered 
magnetization and H e is the so-called exchange field. As 
follows from ©, in this field m =1, i.e. the magneti- 
zation vectors align along the field direction, which is a 
spin-flip transition. After substitution of © and omit- 
ting gradients of m the energy density w in Eq. can 
be written as a function of the vector n alone 



1 

2J l 



[H 2 



(H-n) 2 ]+e Q (n)M 2 



(4) 



For the collinear antiferromagnets, further simplifica- 
tions are possible by restricting the spatial orientation 
of the magnetization vectors. In most practically impor- 
tant cases, the equilibrium configurations m.; and, cor- 
respondingly, the vectors m and 1 remain in or close to 
the plane spanned by the easy axis and the magnetic 
fields In this paper, we chose coordinates with xQz as 
this plane. The restriction on the magnetic configura- 
tions always holds true when the vector H remains in 
the plane spanned by the "easy" and the "intermediate" 
axes of an orthorhombic antiferromagnet. In the uniaxial 
antiferromagnets, in-plane components of the magnetic 
field usually suppress orientational effects of the weak 
in-plane anisotropy . Then, the vectors m, of equi- 
librium states remain in the xOz plane. This means that 
only the uniaxial anisotropy (0) may essentially influence 
the magnetic states. For the magnetic en erg y in terms of 
m and 1, a systematic analysis (see Ref. 0|) shows that 
only the following terms from the uniaxial anisotropy @ 
must be retained, 



,(») _ _ 



= -{K + K')m 2 z - B x l 2 z - B 2 l 4 z 



(5) 



where B\ = K — K', B 2 = (K 20 - K 2l + K 22 )/2. Then, 
within the xOz plane the magnetic states are described by 
only one internal variable, the angle 9 between the easy 
axis and the staggered magnetization. The expansion 
of the energy in terms of 9 with respect to the small 
parameter \e^\/J yields the leading contribution 



w, 



(6) 



(A J)- 1 [{Hi -Hi- H 2 ) cos 26 + 2H X H Z sin 29] , 



H - V2WV/ , 
and the contributions to next order 
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In the transformation of the energy density (0J into the 



simplified energy density wq(9) = Wq{6) +Wq Z> (9) with 
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FIG. 2: (H x ,Hz) phase diagram for easy-axis antiferromag- 
nets. The dashed line H e gives the continuous transition 
into the saturated states. The shaded area is the region 
of the metastable states in the vicinity of the critical field 
(H x = 0,H Z — Ho). Details of the phases diagram in these 
regions are shown in the next figure. 



the contributions from Eqs. © and JSJ), wc have omit- 
ted higher-order terms and a constant combination of 
materials constants. The coefficients in in Eq. © 
are proportional to H 2 . They are generally much larger 



than those in 



(2) 



In wide regions of the magnetic 



phase-diagram 3> and the energy contribution 

(2) 

w Q can be neglected. However, as the magnetic field 
approaches the critical point (H x , H z ) — (0, Ho) the lead- 
ing energy term in © vanishes and the term w 2 ^ from 
Eq. © must be taken into account. In next sections we 
analyze the magnetic states for model J§J, (JHJ. 



III. MAGNETIC PHASE DIAGRAMS 

A. Equilibrium magnetic configurations 

Outside the spin-flop region the equilibrium states are 
described by the energy Wq(9) 10. In easy-axis antifer- 
romagnets the staggered vector 1 interacts with the ap- 
plied magnetic field via the coupling to the magnetization 
vector m _L 1 according to Eq. ©. Correspondingly the 
interaction with the applied field favours the state 1 _L H. 
Note that in isotropic systems one has Ho — 0, and the 
energy density (9) = H 2 cos(2# — 2tp) with ip the an- 
gle between H and z-axis has minima for 9 — if) ± 7r/2. 
In easy-axis antiferromagnets the uniaxial anisotropy ori- 
entates 1 along the easy axis. The competition between 
these two magnetic couplings determines the equilibrium 
states. These states are sketched in the (H X ,H Z ) phase 
diagram, Fig. [3 Minimization of © leads to the well- 



known Neel formula^ for the magnetic configurations 

2H Z H X 



tan 29 = 



m 



H 2 . 



H 2 



(9) 



Because the energy is invariant under the transformation 
9 — » 9 + tt, Eq. © describes solutions with antiparal- 
lel directions of the staggered vector. The equilibrium 
states 9\ t 2 correspond to the minima of the leading en- 



ergy Eq. JHJ), which are given 



(1), 







min[?«Q (9)]. These wells in the potential energy Eq. © 
are separated by a barrier. The height of this potential 
barrier is Aw = max[w(9)] — mhx[w{6)]. For the potential 
© we obtain 



1 

2J 



{HI -HI -H 2 



AHim . 



(10) 



In a magnetic field directed along the easy axis H x = 
the collinear antifcrromagnetic phase for H < Ho and 
the spin-flop phase in the field range Ho < H < H e cor- 
respond to the two branches of solutions for Eq. © with 
9 = and 9 = ir/2, respectively (Fig. ^a),(b)). In mag- 
netic fields that deviate from the easy axis, angular or 
canted phases are realized (Fig. ^c)). The equilibrium 
states described by Eq. © result from the competition 
between the uniaxial anisotropy, which favours the ori- 
entation of the staggered magnetization 1 along the easy 
axis and the magnetic field, which orientates 1 perpendic- 
ular to its direction. At low fields the net magnetization 
is very small |m| <C |1|, so the energy contribution B\l 2 
plays the dominating role for the orientation of the mag- 
netic configuration. The characteristic field Ho given by 
a geometrical mean of the intra-sublattice exchange J 
and the second-order anisotropy for the staggered vec- 
tor B\ measures the scale of the energy contributions 
favouring the easy-axis ground state. Thus, for small 
fields |H| < Ho the anisotropy prevails and stabilizes 
states with the staggered vector 1 nearly parallel to the 
easy-axis, 9 -C 1, with small magnetization, m<l. 

In an increasing field 1 rotates towards the direction 
perpendicular to H, and the net magnetization m grad- 
ually increases. In the (H X ,H Z ) phase diagram (Fig. |2J) 
the hyperbola, H 2 — H 2 = Hq , separates the regions with 
the angles \9\ larger and smaller ir/4. At H = H e the an- 
tiferromagnet transforms into the saturated state with 
rn = 1 along the applied field, m||H. 

The Neel equation Eq. @ has a singularity at 
(H X ,H Z ) = (Q,Ho). At this point (0, Hq) in the phase- 
diagram the main competing forces completely compen- 
sate each other: the leading energy contribution © 

(2) 

equals zero, and the next order energy contribution Wg 
in Eq. © plays a decisive role. In this region, the full 
energy of the model w (0 \9j) can be written as 
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4J 



-Hi cos 2 29 
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FIG. 3: Details of the phase diagram Fig.|5|in the vicinity of the spin-flop field depend on the sign of x as defined in Eq. If ft : 
(a) >«r > 0, (b) yt < 0. Thin lines confine the astroid regions with metastable states. Thick line segments give the first-order 
transition lines. The arrows show spin-configurations in the competing phases along the transition. Hollow points are the end 
points of the first-order transitions. 



with 



Hsf 
AH 



B 2 K 
2Bi 2J' 
H + AH , 



K + K' 
4J 



B 2 
2Bx 



(12) 
(13) 
(14) 



Eq. pi(l represents a consistent expression for the phe- 
nomenological energy near the SF field (0, Ho). It 
includes those higher-order interaction terms that are 
mandatory owing to the compensation of leading energy 
contributions (see Eq. JJJJ). These additional terms con- 
sist of inter- and intra-sublattice uniaxial second-order 
anisotropy terms with parameters K, K' and the fourth- 
order anistropy B 2 of the staggered vector. Comparing 
energy densities Wq 1 ^ from Eq. JHl with Eq. (|ll|l we find 
that these magnetic coupling terms (i) shift the value of 
the SF field H -> H SF = H + AH , where AH < H , 
and (ii) create an additional interaction term propor- 
tional to cos 2 29. This energy contribution stabilizes the 
magnetic states at the SF field (0, Hgp) because it pro- 
vides a potential barrier Awq{Hsf) — \ x \Hq/{2J). 

Due to the relations \K\, \K'\, \K 2k \\ < J, \B 2 \ <C \B X \ 
the region, where and have the same order, is 
restricted to the close vicinity of the point (H x = 0, H z = 
H ): \H Z - H \ < H , \H X \ <C H . In this region the 



energy density wq(0) 



w, 



(i) 



(0) 



- w, 



(2) 



(9) can be reduced 



to the potential expression 

2 J wq sign(j<) 



Ho 



H r cos 2 29 



(15) 



(H z - H SF ) cos 29 + H x sin 29 . 



Here, the characteristic field 



H c 



MHn = 



B 2 
2B 1 



K 
2J 



y / 2JB 1 M < Ho (16) 



sets the scale of the field region around the critical field 
(0,H ), where the interactions described by the energy 
density ifH))) produce noticeable effects. The energy den- 
sity $(0) (|15fl functionally coincides with that of a uni- 
axial ferromagnet 



*/(0) = -^3M o cos 2 



H z cos (j) — H x sin <j> , (17) 



where (3 is the anisotropy constant, <f> is the angle be- 
tween the magnetization M and the Oz axis4£ In the fer- 
romagnetic model (|17fl the anisotropy field H a = \(3\Mo 
corresponds to the characteristic field H c in l)15|l. the 
components of the magnetic field (H x , H z ) correspond to 
(—H x , —(H z — Hsf), and the angle cj> to the angle 29. 
Thus, the behaviour of the uniaxial antiferromagnet is 
reduced to the well-known model and the corresponding 
mathematical results for the magnetic states of uniaxial 
ferromagnets 
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For the model (|15fl the equation for the equilibrium 
state d$/d9 = is 

sign(x).ff c sin 20 cos 29 (18) 
-(H z - H SF ) sin 29 + H x cos 26* = 

and the existence region for metastable states is bounded 
by the astroid 

H 2 J 3 + {H z - H SF ) 2/3 = H 2 J 3 . (19) 

The polar angle for the cusps on the sides of the astroid 

ipc = arctan(|iJ c |/iJ ) = arctan(xr) w x, (20) 

is the so-called critical angle of the spin-flop. This is 
the maximal angle for the existence of the metastable 
states in obliquely applied magnetic fields. The notion 
of this "critical angle" was introduced in Ref. , where 
the astroids of type (|T§)| have been calculated for the 
model with second-order anisotropy. The character of 
the magnetic states within the astroid i|19|) depends on 
the sign of the parameter n <|11|) : 

(i) For k > in magnetic field H z = Hsf, H x = 0, the 
first-order transition occurs between the antiferromag- 
netic and spin-flop phase. This is the proper jump-like 
spin-flop transition. Note that the characteristic field Hq 
in the Neel equation © differs from the spin-flop field 
Hsf, Eq. (|14[) . At finite transversal components of the 
magnetic field H x the first-order transition happens be- 
tween distorted AF and SF, i.e., canted phases at the line 
{H z = H 0} \H X \ < x) (Fig. OJa)). The solutions for these 
competing phases are (cf. Ref.[0|) 

9i = arcsm(H x /H c ), 9 2 = -it/2 - 9 Y . (21) 

For increasing H x < H c the difference between the com- 
peting canted phases decreases. This difference disap- 
pears at endpoints of the first-order transition lines. 
These endpoints are located at the cusps of the as- 
troid dJ), (H c ,H S f), (-H c ,Hsf)- The configuration of 
the magnetic states in these points is for both phases 
0i,2 = T*r/4 (Fig. Eta)). 

(ii) For x < 0, the canted phases exist as stable states 
even in magnetic fields along the easy axis within the 
astroid l|19|l . Minimization of i|ll|) for this case yields the 
deviation angle of the solutions from the easy-axis 

9 b =v-axoaoB[{H z -Ho)/\H D \]. (22) 

These solutions describe a continuous rotation of the 
staggered vector 1 from the AF phase at the low cusp, 
H z = Hq — \H C \, to the SF phase at the high cusp, 
H z = Hq + \H C \ (Fig. Etb)). The plane(s) of this ro- 
tation are determined by in-plane magnetic anisotropy. 
Depending on the crystal symmetry there are several 
such planes (half-planes) spanned by the easy-axis and 
easy magnetization directions in the basal plane. This 
degeneracy of the magnetic states is lifted by a deviation 



of the applied field from the easy axis. In such a field, 
the canted state with the largest projection of m onto 
the field direction corresponds to the stable states, while 
other states preserve metastability for small deviations of 
the field and become unstable for larger deviations. This 
means that the vertical axis of the astroid l|19|) given by 
Ha — \H C \ < H z < Hq + \H C \ is a first-order transition 
line between several canted phases. In particular, in or- 
thorhombic antiferromagnets there is only one plane of 
rotation (spanned by the easy and intermediate axis) and 
the phase transition occurs between two canted phases 
with opposite rotation sense. Such magnetic states have 
been studied in Ref. [l6[ . 

In next subsections we analyze field dependencies of 
the magnetization and the magnetic susceptibility. 



B. Magnetization and magnetic susceptibility 

For the field dependencies of the net magnetization 
M = Aio(nii + 1112) = 2Mom, one derives from 

M x = -^j[H x (l + cos 26») - H z sin 29} 
M z = M^[H Z (1- cos 29) - H x sin 29] , (23) 

where 9 is the solutions of the equations © and i|18|) 
minimizing the energy. For the calculation of the magne- 
tization vector and susceptibility tensor in the vicinity of 
the SF field, it is convenient to rewrite the components 
of the total magnetization 123fl in the following form 

M x = M SF [- sin 26 + ^(9)], 

M z = M SF [(l-cos20)+// 2 (0)] , (24) 

where 

Ixl 

m{9) = y-[H z sm20 + H x (l + cos26)], 
H c 

/j, 2 (9) = i-i [H z (l -cos 26») -H x sm20] . (25) 
H c 

The magnetization 




characterizes the typical values of the net magnetization 
in the spin-flop region. By virtue of the relation H c 3> Hq 
the magnetization Mgp amounts only to a small fraction 
of the saturation value Mq. 

In the SF region, given by H Xl H z < H c , see Eq. JT5J, 
the functions \i\ and /x 2 are very small, /Xj <C 1, i = 1,2, 
and can be omitted. By substituting the solutions l|21|) 
into ll'-il) we obtain the magnetization on the transition 
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FIG. 4: Field dependencies of the magnetization components near the SF field for x > 0: transverse M X (H X , H z )(&), longitu- 
dinal M Z (H X , H z ) (b) and their projections onto the (H z , M x i z \) planes (c), (d). Only the magnetization curves corresponding 
to the stable states are plotted. 



line for x > 0, 



axis, while the parallel components are equal 




Mi 1 ^ = 



(27) 



Mi 1 



M« = M« 




Msf 



H z - H SF 



H z — Hsf 



(28) 



The transverse components M x are equal in both phases. 
The longitudinal components undergo a jump AM Z = 
M SFy /l - (H x /H c ) 2 at the SF (Fig. 0J. The parame- 
ter Msf Eq. (|26|l is equal to the maximum value of the 
magnetization jump at the SF transition. 

For yc < the magnetization in the competing phases 
have antiparallel components perpendicular to the easy 



In this case the transverse components M x undergo a 
jump given by AM, = M SF ^1 - {H z - H SF f / H% at 
the first-order transition. 

The landscapes of the magnetization "surfaces" 
M X (H X , H z ), M X (H X , H z ) near the first-order transition 
and in the adjacent regions of the phase-space are shown 
in Fig.0] The shape of these surfaces is reflected by pecu- 
liarites of the magnetic susceptibility. The components of 
the tensor of the internal static magnetic susceptibility 
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FIG. 5: Field dependencies of \ X x(H z ) for various values of 
i/a;. This and in two following figures display susceptibility 
components for antiferromagnets with yc > 0. The suscepti- 
bility units are l/(4Jx). Only the branches corresponding to 
the stable states are plotted. 




FIG. 7: Field dependencies of Xzz(H z ) f° r various values 
H x . The arrow-like shape of the functions Xxx(H z ) within the 
transition region (\H X \ < H c ) changes into a bell-like shape 
outside this region. In both cases the functions Xzz(H z ) have 
a maximum at H z — Hsf- These maximum values are plotted 
as function of H x in the inset. 




0.0 0.5 1.0 
{H-HJH 



FIG. 6: Field dependencies of Xxz(H z ). Inset (A) shows 

Xiz(ifi) at the transition line H z = Hsf- In Inset (B), the 

location of the extremal values from Eq. are plotted by 

a solid line and their amplitudes from Eq. 1361 by a dashed 
line. 



Xij = dMi j dHj , are derived from the equations (cf . |49j ) 
1 



X.X. 



4J 
1 

** ~ 47 
1 



[1 + cos 29 + 2H 2 n(6) cos 2 (20 - ip)], 



(i) = — fsin 29 + H 2 Vl{9) sin(2(9 - ip) cos(26» - ip)] 



,(i) 



4J 



[1 - cos 26» + 2H 2 n(9) sin 2 (29 - ip)], (29) 



where = d 2 <5>(cos29)/d(cos29) 2 , $(cos26») is de- 

rived from Eq. I|15|) . These relations, together with 
Eq. i|18fl yield field dependencies of the magnetization 



components and the susceptibility tensor in the spin-flop 
region. Near the SF field the expansion of l|29|) with re- 
spect to the small parameter xr <C 1 allows a considerable 
simplification of the expressions for Xij , 



1 



sin 26»cos 2 29 



Xxj 



4J|x| (sign(x) sin 3 20 + # x /if c ) 

1 cos 3 29 

"Ijjxj (sign(x)cos 3 29 -H z /H c ) ' 



(30) 



1 



sin 2 29 cos 29 



4J|>r| (sign(x) sin 3 20 + /J x /.ff c 
1 sin 20 cos 2 20 



4J|x| (sign(x) cos 3 29 - H z /H c ) ' 



(31) 



1 



sin 3 29 



4J|x| (sign(x) sin 3 29 + H x /H c 
1 sin 2 29 cos 29 



4J|x| (sign(x)cos 3 26»- J ff z /7J c ) 



(32) 



For antiferromagnets with x > the expressions from 
Eqs. i|3U|) . (|3*T)l . and (|32*)l for %y m dependence on H x are 
convenient for numerical calculations and those depend- 
ing on H z in the case x < 0. As the field approaches the 
stability limits, given by the lines $eg = d 2 &/d9 2 — > 0, 
the functions Xij ^ l/^ee diverge. For applications to 
real antiferromagnets it is important to distinguish the 
branches of the functions Xij{H x , H z ) corresponding to 
the stable states. Typical field dependencies of the com- 
ponents of Xij for the stable phases are plotted in Figs. 03 
El and0 The functions Xxx(H z ) have qualitatively differ- 
ent field dependencies within and outside the region near 
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the first-order transitions. In the former case \H X \ < H c , 
the curves Xxx(H z ) have arrow- like shape with a max- 
imum at the transition field H z — H$f- In the latter 
case \H X \ > H c , the curves Xxx(H z ) have a minimum in 
fields along a line prolongating the transition line and 
two symmetric maxima 



(max) 
A-xx 



I 



1 - 



1 



3 cos 2 a 



(33) 



for fields lower and higher than the SF field i?sF- The 
location of these maxima in (H X ,H Z ) for H x > H c is 
determined parametrically by a set of equations 



H 



(max 1) 



H 



(max 1) 



2H C cos a 
3 cos 2 a — 1 
H c sin 3 a 



3 cos^ a — 1 

with a = —tt/2 — 29. For \ xz (Fig. stationary points 

1 3cos 2 a-2 

(35) 



ixir x) i = 



8Jx 31 sin a I cos a 



are situated on the lines 



j_j-(m&x 2) 



H c cos 3 a 



Hi 1 



= ± 



3 cos 2 a — 2 ' 
2H C sin 3 a 



3 cos 2 a — 2 



(36) 



Similarly, functions dxij/dHk(H) also show qualitatively 
different dependencies in different parts of the phase di- 
agram and include extremal points. The field values for 
all these anomalies are determined by the material con- 
stants of the antiferromagnets. This connection between 
anomalous field-dependence of the susceptibility tensor 
in the SF region and materials paramters can provide 
a basis for experimental approaches to investigate the 
magnetic properties and materials parameters in antifer- 
romagnets. Values of magnetic interactions, in particu- 
lar higher-order contributions to the magnetic anisotropy, 
can be determined by measurements of the susceptibil- 
ity tensor in the SF region and fitting of the Xij data 
to the theoretical expressions shown above. This will be 
demonstrated in section Ivl 

Here, we compare characteristic values of the mag- 
netic susceptibility for antiferromagnets in the major part 
of the magnetic field phase diagram and in the SF re- 
gion. In the (H X ,H Z ) phase diagram (Fig. 0) saturation 
is achieved at the exchange field \H\ = H c . Thus, the 
average susceptibility is (x)afm ~ Mq/H b = 1/(2 J). 
This is exactly the value of the susceptibility in the SF 
phase. The metastable region in the phase diagram near 
the SF field has a width AH = 2H C = 2H i/) c , Eq. JSJ, 
while the magnetization changes by AM = MqHq/ 'H e> 
Eqs. (|27[l . (|28[l . Thus, the average susceptibility in this 
region can be estimated to (x)sf = M)/(2|x|i? e ) = 
1/(4 Jh) ~ 1/(4 Jipc)- This average equals the expres- 
sions from Eqs. I|3U|) . I|31|) . and (|32|1 up to some numeri- 
cal factor. The magnetic susceptibility near the SF field 



is strongly enhanced compared to the average suscepti- 
bility in the major part of the phase diagram, (x)sf = 
(x) afm (Ho/(2H c )). This enhancement is given by the 
ratio (Hq/(2H c )) = l/(2ip c ) which usually amounts to 
a factor of several hundreds. Note that the absolute 
change of the magnetization AM = MqHq/H c <C Mq 
is tiny. However, due to the extremely narrow width of 
the metastable region near the SF field, the magnetic 
susceptibility becomes very strong in this region. 

As was mentioned above, similar functional expressions 
describe the field dependence of the magnetization M^) 
for the model of a uniaxial ferromagnet (|17|) and for the 
antifcrromagnet in the SF regions following the basic ap- 
proach Eq. I|24|l with fa = 0. In fact, one can derive 

M(H/H C ) =(^V (/) (H/ffa) + M (||V (37) 



2H r 



The last term in (|37J) signifies a shift in direction of the 
easy axis a. A corresponding relation between the com- 
ponents of the magnetic susceptibility near the SF tran- 

( f) 

sition and those of the ferromagnetic susceptibility x\j 
is given by 



Xij(K/H c 



H, 



2H r 



x\hn/H a ). 



(38) 



These relations demonstrate the physical similarity of the 
field dependencies for the magnetic properties in uniax- 
ial ferromagnets and in easy-axis antiferromagnets near 
the SF field. This equivalence is established by introduc- 
ing reduced units and a shift for the magnetization in 
Eq. 1(57} . and by reduced units and a scale- factor for the 
susceptibility in Eq. I|38|l . It is a consequence of the for- 
mal similarity in the phenomenological models Eqs. 1|17|) 
and l|24|) for both systems. The relations Eqs. I|37|l and 
(|38|l are useful for comparative studies on reorientation 
transitions in uniaxial ferromagnets and antiferromag- 
nets near the spin-flop. 



C. Two-scale character of easy-axis 
antiferromagnets 

The magnetic states analyzed in the previous two sec- 
tions display the pronounced two-scale character of easy- 
axis antiferromagnets. The magnetic phase diagrams in 
Fig. O comprise the main features of the solutions in the 
spin- flop region. On the large scale of the (H X ,H Z ) phase 
diagram two materials parameters rule the behaviour: 
the exchange or saturation field H e and the SF field 
H defined in Eqs. and J7J), respectively (Fig. 0). 
The former characterizes the strength of the antifcrro- 
magnetic exchange interaction, the latter comprises the 
interactions favouring easy-axis states. The equilibrium 
orientations of the staggered magnetization, derived from 
the Neel Eq. @, result from the competition between 
these interactions and the applied field. Thus, the field 
Hq sets the characteristic scale in the major part of the 
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magnetic phase diagram. Note, this is the only material 
parameter in the Neel equation JJjJ. 

In the vicinity of the SF field, where two of the main 
energy contributions cancel each other, much weaker (rel- 
ativistic) interactions enter the scene set by Eq. Ql lfl . The 
characteristic scale in this region is given by the field H c , 
as introduced in Eq. i|16|) . Hence, \H C \ <§C Hq defines a 
"fine" scale of the system. This scale gives the value of 
the potential barrier Awo(Hsf) — H C M at the spin-flop 
field and it fixes the size of the metastable regions around 
the field H , i.e. the astroids. The field \H C \ includes 
two physically different contributions. One of them is 
the ratio between the sublattice second-order anisotropy 
K and the antiferromagnetic coupling. The other is the 
ratio between fourth- and second-order anisotropies, Bi 
and Bi, of the staggered vector (see Eq. |JSJ). Generally, 
the two terms have same order of magnitude. 

The strengths of the magnetic couplings in usual 
antiferromagnetic materials obey a well-defined hier- 
archy with very strong exchange and weaker uniaxial 
anisotropy. This hierarchy is given by the relations J 3> 
K,K',B 1 » B 2 . Hence, the field H = H cy /B 1 /(2J) 
is much smaller than the exchange field H c , and H c = 
(K/J + B2I Bi)Hq/2 is much smaller than Hq. Corre- 
spondingly, the jump of the magnetization at the SF tran- 
sition is small Am ~ H /H o < 1 (see Eqs. pTjl. (j2Hl)- 
The potential barrier separating the stable states at the 
SF field Aw(Hsf) = kB\Mq is again much smaller than 
the barrier in the ground state, Au>(0) = B\Mq, and 
the region of the metastable states is restricted to a close 
vicinity of the SF field, AH ~ H c . This causes the un- 
usually high sensitivity of the magnetic states near SF 
field with respect to small changes in strength and direc- 
tion of applied fields. For example, in an antiferromagnet 
with k > the rotation of the magnetic field H = Hsf 
from ijj = —ip c to ip — ip c causes a change of the stag- 
gered magnetization from 9 — it/ '4 to 8 = — ir/4. This 
lability of the magnetic states is the underlying reason 
for the remarkable magnetic effects in the SF region. 



fields breaks down at field-induced phase transitions>i2i^ 
In the vicinity of such transitions the homogeneous states 
are unstable with respect to the transformation into mul- 
tidomain structures consisting of domains formed from 
the competing phases^S^ Within the thermodynamic 
approximation or generalized phase theory, two-phase 
multidomain states are described by the equation 



H (o) = H tr + 4ttN < M >, 



(40) 



where H tr is the field value for the first-order 
transition.!^ Here, < M >= frM^ + ^M' 2 ', M; 
is the total average magnetization, and the variable pa- 
rameters £j are the volume fractions of the coexisting 
phases (i — 1,2, £1 +£2 = 1). For the SF transition the 
phase theory approximation has been found to be valid 
practically in all regions of the phase diagram, where 
multidomain states exist>i2ii& 

For model (|15f) with xr > the phase theory equations 
(|40|l can be written in the following form 

= H x (l + r)N xx ) (41) 
+r)N xz (& - + V N XZ H C , 

= Hgp + H x fjN xz 

+vN zz (£i - &)y/H*-H2 + V N ZZ H C , 

where H x varies along the first-order transition line 
( \H X \ < H c , H z = Hsf)- The parameter 77 = 
2ttMsf/H c — n/(>cJ) measures the ratio of the stray field 
energy and the potential barrier between the coexisting 
phases at the SF transition. Equations (|41|> allow us to 
derive the parameters of the multidomain states, H x , £j, 
as functions of the external field. In particular, for the 
relevant case with N xz = 0, Eas. H41|) with £]/ 2 ) = 1 yield 
the boundary of the multidomain states as an ellipse 



H 



(e) 



H c + 2irN xx M SF 



H 



(e) 



H SF - 2irN zz M SF 



2irN zz M SF 



(42) 



IV. DEMAGNETIZATION EFFECTS AND 
MULTIDOMAIN STATES 

In the previous section the equilibrium magnetic con- 
figurations have been derived as functions of the internal 
magnetic field H, which differs from the applied external 
field H^ 6 ) due to the demagnetization field H^ m ' of the 
sampled In a homogeneneously magnetized ellipsoidal 
sample (including the limiting cases, i.e. plates and long 
cylinders) the equation 



H (o) = H + 4ttNM(H), 



(39) 



establishes the relation between the external and inter- 
nal magnetic fields. This relation allows to express the 
solutions for magnetic states as functions of the external 
field by using the demagnetization tensor N of the sam- 
ple. However, the relation between internal and external 



with semiaxes a = H c + 2itN xx Msf, b = 2itN zz Msf 
(Fig-IBIA)). For x. < similar equations yield the bound- 
ary ellipse 



r h-w " 
n x 


2 r 

+ 


2ttN xx M sf 





(e) 



H SF - 2irN zz M SF 



(43) 



H c + 2irN zz M SF 

with semiaxes a = 2itN xx Msf, b = H c + 2itN zz Msf 
(Fig. Hp)). 

The largest tilt angles between the applied field and 
the easy axis, at which the multidomain states still exist, 
can be readily derived from Eqs. I|42(l . I|43(l . 



. (o) H c + 2irN xx M SF 

%i = 



H 



ttN x 



SF 



,( c ) _ 2kN xx M S f 
^ u ~ H SF 



kN x 
J 



J 

for >c < . 



for k > . 



(44) 
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FIG. 8: Magnetic phase diagrams in components of the external field (Hi ; , -ff* ) for x > (a) and x < (b) include regions 
of the multidomain states (shaded areas) . 



The phase diagrams in Fig. Q3] demonstrate the strong 
formal resemblance for the two qualitatively different 
cases. In both cases thermodynamically stable multido- 
main states arise near the SF transition in a close vicinity 
to the SF field. However, due to the different character of 
the phase transitions for x > and x < 0, Eqs. (|21|l and 
Effi . the evolut ion of the magnetic states within the mul- 
tidomain regions H42|) and (|43[) is different. For h > 
the competing phases <|21|) coexist along the horizontal 
line segment (H z = H SF , \H X \ < H c ). In the (H^ , il£ e) ) 
phase plane a set of the vertical straight lines describes 
regions with these fixed transition fields (Fig. [Hj(a)). Due 

to the smallness of ipty 144|) external magnetic fields with 

fixed directions (|V ;< ' C ' ) | < ipci practically intersect the re- 
gion of the multidomain states l|42|) along the lines with 
a fixed transition field. A variation of fields with such a 
fixed orientation causes in the system magnetization pro- 
cesses through the displacement of domain walls between 
the coexisting states. On the other hand, rotating mag- 
netic fields with fixed amplitude that cross the region i|42|) 
mainly cause a continuous deformation of the magnetic 
configurations within the domains. For x < rotating 
fields cross the multidomain region 143(1 almost perfectly 
along lines corresponding to fixed internal fields. Thus, 
a rotating field induces a redistribution of the domains. 
On the other hand, magnetic fields with fixed directions 
produce mainly reorientation effects within the domains. 



Correspondingly, the limiting angles 144fl have a dif- 
ferent physical meaning in both cases. For systems with 
k > the critical points at ip^ = i^ci^i = 



Hsf + 2ttN zz Msf correspond to internal states where the 
difference between the magnetic configurations in the co- 
existing phases disappears. For antiferromagnets with 
k < the transition into the homogeneous state at the 
critical points t/A) = i^n, # (c) = H SF + 2irN zz M SF oc- 
curs by a complete replacement of one of the coexisting 
phases by the other through domain processes. 

The quantitative description of the magnetization in 
the multidomain states is provided by Eq. (|4U|I . This 
equation together with those for M(H) lt2"71> . (|2"5t yields 
the functions (M)(H( C )). It follows from gOJl that the 
magnetization is a linear functions of the applied field 
only when H^ ) varies along lines corresponding to fixed 
internal transition fields. In the general case the variation 
of the transition field HW(H( C ') causes complex depen- 
dencies of the magnetization on the external field in the 
multidomain state. Differentiation of (M) in Eq. 14U|) 
yields the equation for the components of the magnetic 
susceptibility (cf. 0) 



47T7V, 



ijXjk 



= 6. 



(t) 



ik 



dH, 



(45) 



The first term on the right side of l)45|l describes the 
process of a redistribution of volume fractions for the 
different phases through displacement of domain walls, 
while the second term is associated with the variation of 
the magnetic states within the domains by rotations of 
the magnetization. 

In homogeneous phases Eq. (|39fl allows one to relate 
the external susceptibility arid the internal suscepti- 
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FIG. 9: Dependencies of xi°2 on external field calculated for (^HsNHs^CuCU for different fixed values of ip^: 
0.1°(i),1.50 o (2),2.86 o (3),2.92 o (4),2.98°(5),3.50°(6). Black solid lines correspond to the homogeneous phases and dashed 
(blue) lines to the multidomain states. Hollow points indicate the boundaries between these regions. Variations of the internal 
magnetic field H(H^) for fixed directions of the external field, ip^ = const (b). The (red) arrow in (b) indicates the location 
of the endpoint of first-order transition. The curved "trajectories" H(H^ explain the reentrant character of the X^x(H^) 
functions in (a). 




FIG. 10: Experimental dependencies of Xxx on the external 
field for a number of fixed angles tp^ for a sphere of 

(C 2 H 5 NH3)2 CuCU (Sample No. 1) at T = 4.2 K. 

bility x 

X (o) - (V 1 +4vrJV) 1 . (46) 

Eqs. I|4(J|) . H45fl . <|46[l transform the components of the 
net magnetization (|24() and magnetic susceptibility (1301) . 
ijjnj ! H^2(l derived in the components of the internal field 
into those in components of the external field. The 
tensor of external susceptibility l|4bl includes the in- 
ternal susceptibility x an d the stray field contribution 
Xm = (47rA r )~ 1 known as shape susceptibility. 10 In the 
multidomain state with fields varying along the lines of 
the fixed transition field the external susceptibility as 



given by l|45|l has only contributions from the shape sus- 
ceptibility, 4:TrNijXjk — Sik- When the evolution of a 
multidomain state involves a variation of the transition 
field a specific susceptibility contribution arises that is as- 
sociated with the rotation of the magnetic states within 
in the domains. In particular, for the multidomain states 
in Fig. |S] along the line Hx = 

X« = (*>0), (47) 

For x > the magnetic field varies along the line of the 
fixed transition field, and the susceptibility includes only 
shape contribution. For x < the magnetic field varies 
along the line with fixed volume fractions, £i = £2 = 1/2, 
and the evolution of the system consists of a continuous 
reorientation in the domains i|22|) . For this process the 
net magnetization is derived from Eq. 1)28(1 . and the in- 
ternal susceptibility is Xzz = V(4|>f|J). The external 
susceptibility Xzz i|47|l includes both internal and shape 
contributions. Generally during the SF transition the 
values of the internal susceptibility l/(\x\J) H30|) . (13 1 1) . 
(|32|l arising due to variation of the homogeneous mag- 
netic states and the shape susceptiblity 1 /N originating 
from the reconstruction of the multidomain states are of 
the same order and larger than the values for the external 
susceptibility outside of the SF region. For example, in 
the external field along the easy axis x*^ equals zero in 
the AF phase, and in the SF phase xfl = 1/{2J+AttN zz ). 

Fig. El (a) shows the calculated external-field de- 
pendencies of Xxx{H^ e ') for a number of fixed angles 
^M e ) and with materials parameters close to those for 
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FIG. 11: Location of extremal points of Xxx and dxxx/dH^ 
in the phase plane (H (c) , t/> (o) ) for an elliptical (C2H 5 NH 3 )2 
CuCLt sample (No. 3). The shaded area indicates the region 
of multidomain states. 



(C2H5NH3)2CuCi4. According to previous investiga- 
tions in this antiferromagnet, xr > 0. 20-49 The functions 
Xxx(H^) for small angles ip^ < ipty , are given by lines 
1-3 in Fig. 0(a). They cross field-dependent regions cor- 
responding to homogeneous phases (black solid lines) and 
areas with constant susceptibilities (dashed blue lines) 
indicating multidomain regions. For V (e) < 4>a the 
functions Xxx(H^) (lines 5, 6) have characteristic fea- 
tures as described in the previous section (Figs. [S] and 
Eqs. (|33|l . (|34|) h The intermediate line 4 shows ambiva- 
lent features. In the center it includes a jump into the 
multidomain states as the lines for ip^ < , but out- 
side this narrow central feature the function \xx 
behaves similar to those for ip^ > ipj$. This interesting 
effect is explained by strong deviations of the total mag- 
netization from the easy axis in the vicinity of SF field. 
The enhanced values of the tranverse magnetization M x 
(Fig. 0| (c)) create a strong demagnetization screening 
near the SF field. Accordingly, for the external field vary- 
ing along lines ip^ e ' = const, interal field "trajectories" 
H(H^) deviate towards the H z axis near the SF field, 
as depicted in Fig.|§l(b). The trajectory for the evolution 
of the internal field H starts from outside the transition 
region i.e. \H X \ > H c but enters this area in the vicin- 
ity of the SF field for the external field applied under an 
angle ^ (c) = 2°. 



V. MAGNETIC PHASE DIAGRAM OF 

(C 2 H 5 NH 3 )2CUCL4 

An experimental investigations of the differen- 
tial magnetic susceptibility was carried out for 
(C2H5NH3)2CuCl44£ This model antiferromagnet has 
orthorhombic lattice structure with space group is Pbca, 
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V (e) (deg) 

FIG. 12: Experimental phase diagram in components of ap- 
plied field for (C2H5NH3)2 CuCU (for spherical and ellip- 
soidal samples). The region of multidomain states is within 
the area circumscribed by measured points. Filled circles give 
the field values for the maxima in d X xz/dHi c) (Hi c) ); hollow 
circles give the inflection points in dxzz/dH l f\lli e ' > ). 



and orders at a Neel temperature Tn — 10.20 K£. The 
magnetic structure consists of ferromagnetic layers par- 
allel to the xOz planes in the notation of this paper 
with weak antiferromagnetic couplings. According to 
Ref. Q the ferromagnetic intra-layer interactions corre- 
sponds to an effective field 500 kOe, while antiferromag- 
netic coupling between layers is JM = 837.5 Oe. A 
weak second-order anisotropy in the xOz plane KM = 
76 Oe stabilizes the collinear ground state with the sub- 
lattice magnetizations along the z-axis. At T = 4.2 
K the spin-flop field is at Hgp = 305 Oe. Magnetic inter- 
actions in (C2H5NH3)2CuCi4 include a Dzyaloshinsky- 
Moriya coupling described by an energy contribution 
wd = D(M\ x M2 y — Mi x M\ y ) with an effective field 
DMo= 119 Oe£ Generally this interaction rotates the 
magnetization vectors away from the basal plane xOz. 
However, due to the strong orthorhombic anisotropy with 
a value of 1504 Oe, the vectors M.; are practically con- 
fined to the basal plane. Deviations from this plane do 
not exceed 2°i2i Thus, the planar model l|ll|l can be ap- 
plied to this antiferromagnet. 

The differential magnetic susceptibility components 
were measured with an inductive technique using three 
pairs of modulating and pick-up coils arranged along per- 
pendicular directions. The modulating fields had ampli- 
tudes between 0.3 and 1.0 Oe at frequencies of 9, 17- 
21, 86 and 133 Hz. The three samples used for the 
measurements were cut from a (C2H5NH3)2CuCU sin- 
gle crystal. The samples had the following geometrical 
parameters: (Sample No. 1) a sphere with diameter 3.1 
mm (N xx = N zz = 1/3); (Sample No. 2) an ellipsoid 
with axes a x = 5.00 mm, a y — 1.75 mm, a z — 5.00 mm 
(N xx =0.185, N zz = 0.185); (Sample No. 3) an ellip- 
soid with a x = 3.43 mm, a y = 2.45 mm, a z — 5.75 mm 
(N xx =0.344, N zz = 0.186). All the measurements of 
components Xij have been carried out at T = 4.2 K. 
The recorded data for all investigated samples are in 
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close accordance with the theoretical results expounded 
in previous sections. An example is shown in Fig. (|1U[). 

where \xx components are plotted for the spherical sam- 
ple. The experimental data follow closely the theoretical 
results of Eq. lIHUjl as sketched in Fig. |5| In particular, 

the dependencies of (# (c) ) for V> (c) = 3.5°, 4.0° and 5° 
display the reentrant behaviour imposed by the rotation 
of the internal field towards the ff z -axis, as demonstrated 
in Fig. |SJ By fitting experimental data for the com- 
ponents and their field derivatives with the theoretical 
dependencies, values of the material parameters can be 
deduced and the magnetic phase diagram of this antifcr- 
romagnet has been constructed. Fig. 1111 shows the loca- 
tions of the extremal points for xifi and d\xx /dH^ and 
the region of the multidomain states in the (H^ ,ip^) 
phase diagram of the ellipsoidal samples. 

The experimental (H^^ip^) phase diagram for the 
spherical sample is shown in Fig. ^| In particular, for 
the maximal angle -0^' the following results have been 
obtained (spherical sample No. 1) 3.50° ± 0.65° (sample 
No. 2) 2.70° ± 0.45° (sample No. 3) 3.55° ± 0.65° These 
results yield the values of the spin-flop field Hsf = 306 
Oe, the critical angle ip c — 1.7°, and the characteristic 
fields H = 300 Oe and H c = 9.1 Oe. 

VI. CONCLUSIONS 

Magnetic configurations for an extended class of two- 
sublattice collinear easy-axis antiferromagnets have been 
obtained as functions of the values and directions of 
the applied magnetic field, and the corresponding mag- 
netic phase diagrams have been constructed. The mag- 
netic behaviour of these materials strongly depends on 
the strengths of the applied field. In the main part of 
the magnetic phase diagram they are described by the 
well-known Neel model (JfoJ). In the vicinity of the SF 
field they are described by an effective model for a re- 
orientation transition that is formally equivalent to the 
Stoner-Wohlfarth model of a uniaxial ferromagnet (|15fl . 



This two-scale character of easy-axis antiferromagnets 
has been ignored in most previous investigations. The 
analysis of magnetic-field-driven reorientation effects and 
the concomitant multidomain states provides a consis- 
tent picture of the magnetization processes near the SF 
transition. At the SF transition the weak intrinsic higher- 
order couplings of the antiferromagnetic material and the 
dipolar stray field cause important and noticeable effects 
and must be included in the analysis. These effects are 
responsible for the prominent anomalies in the magne- 
tization curves and the magnetic susceptibilities of anti- 
ferromagnets in the spin-flop region. The anomalies pro- 
vide effective approaches for investigations on magnetic 
properties of antiferromagnets that are usually hidden 
away from the spin-flop region. The validity of this ap- 
proach has been demonstrated by an application to an 
orthorhombic antiferromagnet. 

The results on bulk antiferromagnets may also be ex- 
tended to confined antiferromagnetic systems by includ- 
ing surface/interface-induced interactions into the phe- 
nomenological models. In Ref. [5(j it was shown that the 
interplay between surface-induced and intrinsic magnetic 
interactions yield a rich variety of specific magnetic states 
including spatially inhomogeneous twisted states in the 
vicinity of the SF field. The further development the the- 
ory for such cases will be important for an understanding 
of the magnetization processes in ferro/antiferromagnetic 
bilayerS) 3 ^^^ and in antiferromagnetic nanoparticles^ 

The phenomenological model of the two-sublattice an- 
tiferromagnet and its variants can be adopted 
also to describe magnetic states in synthetic 
antiferromagnetsi^i 3 ^ 
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